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Abstract 

We investigate topological charge and the index theorem on finite lattices 
numerically. Using mean field improved gauge field configurations we calculate 
the topological charge Q using the gluon field definition with C(a 4 )-improved 
cooling and an 0(a 4 )-improved field strength tensor F^ u . We also calculate 
the index of the massless overlap fermion operator by directly measuring 
the differences of the numbers of zero modes with left- and right-handed 
chiralities. For sufficiently smooth field configurations we find that the gluon 
field definition of the topological charge is integer to better than 1% and 
furthermore that this agrees with the index of the overlap Dirac operator, i.e., 
the Atiyah-Singer index theorem is satisfied. This establishes a benchmark 
for reliability when calculating lattice quantities which are very sensitive to 
topology. 
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I. INTRODUCTION 



The connection between the topology of a background gauge field and fermion zero modes 
is related to the axial anomaly and the large r\ — t( mass splitting in QCD. Lattice gauge 
theory is the tool best suited to the study of these nonperturbative issues. We study this 
connection in QCD formulated on a periodic lattice, i.e., on a four-dimensional toroidal 
mesh. For a fine enough lattice and/or a sufficiently smooth gauge field configuration we 
should recover the results for continuum QCD on a 4-torus. In particular, we should recover 
the Atiyah-Singer index theorem Jj]]. 

In the continuum the Dirac operator, D = / -f fJ- (d tl +igA tl ), of massless fermions in a smooth 
background gauge field with nontrivial topology has eigenmodes with zero eigenvalue (i.e., 
"zero modes" ) which are also chiral eigenstates of positive or negative chirality. The Atiyah- 
Singer index theorem |1J gives the result 

Q = index(X>) , (1.1) 

where 

Q = 32^/ ^WM^^) (1-2) 

is the topological charge of the background gauge field and where 

index(.D) = n_ — n + (1-3) 

is the chirality index of the Dirac operator. Here n + and n_ are the number of zero eigen- 
modes with positive (right-handed) and negative (left-handed) chiralities respectively, i.e., 
Dijj = with +/— chiralities such that 75^ = =b0. These results apply to QCD defined on 
a continuum 4-torus, where ultimately we wish to take the size of the 4-torus to infinity, i.e., 
the infinite volume limit. 

However, in the lattice formulation one has ambiguities associated with the discretiza- 
tion and in general one can only expect the index theorem of Eq. (|1 . 1| ) to be satisfied on 
a sufficiently fine lattice and/or after sufficient smoothing of the gauge fields. The use of 
improved operators and actions will lead to the index theorem being satisfied for less strin- 
gent conditions on lattice spacing and/or smoothness. There are different ways to calculate 
Q for a given gauge field configuration, e.g., using the gauge-field tensor definition of 
Eq. (|1.2|) or using the geometrical definition In this study we focus on the definition in 
Eq. ( (L.2[ ). On a continuum 4-torus the two definitions of topological charge for the gauge 
field are necessarily identical. However, in the calculation of hadronic observables using 
typical lattices, the configurations in the ensembles are too coarse for lattice definitions of 
Eq. ( |1.2| ) to lead to integer topological charge and so the index theorem is not satisfied. 

In an arbitrary gauge field, the Wilson fermion operator does not have exact zero modes 
due to the Wilson term which removes fermion doublers and breaks chiral symmetry. At- 
tempts to study the index theorem with such an action requires one to estimate the number 
of "zero modes" by looking at low lying real eigenvalues M. Due to the difficulties in esti- 
mating the index of the Wilson fermion operator in a precise manner, it is difficult to reach 
a definite conclusion concerning the validity of the index theorem on a finite lattice with 
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Wilson fermions. In the case of Ginsparg- Wilson (GW) Dirac fermions |]5|-[7| there is an 
exact lattice realization of chiral symmetry || and the GW Dirac operator possesses exact 
zero modes. Hence the ambiguity associated with the need to subjectively estimate the 
number of "zero modes" of the Wilson fermions is absent. A good review and introduction 
to many of these issues can be found in Ref. ||. 

There are several numerical studies |T0|-[T5| of the index theorem on a finite lattice. In 
Ref. |T(| the overlap fermion formalism is used to estimate the probability distribution of 
topological charge p(Q) in pure SU(2) gauge theory by examining the spectral flow of H(/i), 
where H (/x) is the hermitian Wilson-Dirac operator. The study in Ref. f|TT|] explores the real 
eigenvalues of the Wilson-Dirac operator, which are identified as the lattice counterparts of 
the continuum zero-modes. They studied topology by examining the complete spectrum of 
the fermion matrix in SU(2) gauge theory on small lattices. Ref. |]12[ uses the spectral flow 
method to perform a comprehensive study of both quenched SU(3) and dynamical fermion 
configurations. There the role of the SW term is also examined. 

We focus here on Neuberger's overlap Dirac operator, which is an explicit solution of the 
Ginsparg- Wilson relation, and investigate the Atiyah-Singer index theorem on a finite lattice 
by numerical methods. In a recent paper ||16|| , Adams performed an analytical study which 
showed that the index theorem was satisfied by the overlap operator in the continuum limit. 
We calculate the index of the overlap fermion operator directly by measuring the number of 
left minus right zero modes of the overlap Dirac operator D. 

In measuring the gauge-field definition of the topological charge, Q, we use an C(a 4 ) 
improved definition of the field strength tensor leading to an improved topological charge 
operator. We also use a mean-field improved Symanzik gluon action in the quenched ap- 
proximation to generate our ensemble of configurations. Where we have employed cooling 
to smooth our configurations, we have used an 0(a 4 ) improved gluon action in the cooling 
algorithm. The resulting Q approaches integer values after a few cooling sweeps and has 
been verified to be stable for hundreds to thousands of cooling sweeps. 

The paper is organized as follows: In Sec. [H| we introduce the details of our calculation 
of the improved gauge-field topological charge Q. In Sec. |H| we review the overlap quark 
propagator and describe our calculation of index(D). Our results are described in Sec. [TV] 
and finally our summary and discussions are presented in Sec. [V|. 



II. IMPROVED COOLING AND TOPOLOGY 

To investigate the topology of gauge fields on the lattice we first construct an ensemble of 
gauge field configurations. In lattice QCD, the gluon fields are represented by SU (3) matrices 
on each link connecting adjacent lattice sites. These links are parallel transport operators 



between the lattice sites. We use a parallel Cabibbo-Marinari |17| pseudo-heatbath algorithm 



with three diagonal SU(2) subgroups looped over twice and appropriate link partitioning 



For typical lattice spacings used in simulations of QCD the link configurations in the 
ensemble have fluctuations on many scales. In particular, such typical link configurations 
are not smooth at the scale of the lattice spacing. However, the Boltzmann factor e~ Saau3e 
with S gauge 1/g 2 will ensure that link configurations become increasingly smooth as the 
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continuum limit is approached (g — > 0, a — > ). But the cost is that the volumes that one 
can afford to simulate at also become correspondingly smaller. 

One common approach to probing the medium to long range topological structure of typ- 
ical link configurations is to cool them sufficiently that they become approximately smooth 



at the scale of the lattice spacing |19| . Cooling involves recursively modifying the link values 
to locally minimize the action. As we sweep repeatedly over the whole lattice we reduce 
the total action, smoothing out fluctuations on successively larger scales ||20fl . This quickly 
eliminates the high-frequency, rough components of the fields, leaving relatively smooth 
topological structures. Once the link configuration has been cooled sufficiently, we can ex- 
pect the cooling process to preserve the global topological charge Q. If we cool indefinitely 
a link configuration on a sufficiently fine lattice such that we preserve the global topological 
charge Q, then we should converge toward the minimum action solutions with charge Q, 
i.e., we should converge to self-dual configurations (F^ u = F^ u ). The resulting self-dual con- 
figurations on the 4-torus should contain only \Q\ instantons (if Q > 0) or anti-instantons 
(if Q < 0). The action of an instanton or anti-instanton in infinite space-time is known 
analytically and is given by So = 8ir 2 /g 2 . Instantons and anti-instantons have the property 
that they are scale invariant, i.e., they can have any scale in the infinite 4- volume continuum 
and their action So is unaffected. In infinite space-time in the continuum, a self- dual gluon 
configuration with action S will necessarily have topological charge \Q\ = S/S and contain 
only instantons (or anti-instantons). To the extent that the finite volume of a 4-torus is 
sufficiently large compared with the size of the instantons on it, one can expect the above 
results to hold . 

However, it is well established that cooling with the standard Wilson action eventually 
destroys (anti-)instanton configurations, due to the discretization errors inherent in the 
Wilson action pi| . The Wilson gluon (i.e., Yang-Mills) action at each lattice site is calculated 
from the plaquette, a closed product of four link operators 

Swii =/?EE^etr(l- U, v ) (2.1) 

where the plaquette operator is 

= Up(x)U v (x + n)Ul(x + u)Ul{x) . (2.2) 

This Wilson plaquette action contains deviations from the continuum Yang-Mills action of 
O (a 2 ) where a is the lattice spacing. This problem may be remedied by improving the action. 
Tree- level improvement of the classical lattice action combined with mean- field (i.e., tadpole) 
nonperturbative improvement |EJ provide a simple and very effective means of eliminating 
lowest-order discretization errors. The simplest improvement of the classical gauge action on 
the lattice is achieved by taking a linear combination of the Wilson plaquette and an a x 2a 
rectangle, i.e., Symanzik improvement. Since the plaquette and rectangle have different 
0(a 2 ) errors they can be added in a linear combination in such a way that these 0(a 2 ) 
errors cancel leaving only errors of 0(a 4 ) in the classical (i.e., tree-level) gauge action. 
One can preserve this improvement at the nonperturbative level by combining this with 
mean-field improvement of the links ||23|| . More recently DeForcrand et al. [24| have used 



tree-level improvement to construct a lattice action which eliminates 0(a 4 ) errors and leaves 
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only O (a 6 ) errors, by using combinations of up to five different closed- loop products of link 
operators (Wilson loops). From these five loops, particular linear combinations were studied 
in detail and were referred to as 3-loop, 4-loop, and 5-loop improved actions. The difference 
is the number of non-zero contributions in the linear combination of the five planar loops. 
They represent the improved action as 

5 

•Simp = 2_j C i^i i (2-3) 
i=l 

where the Si are the actions calculated as per equation (|2.1j) , but using five different Wilson 
loops, and the q are the weighting constants tuned to eliminate the 0(a 2 ) and 0(a 4 ) errors. 
Additional details may be found in Ref. [f24 . 



In our work we construct improved actions utilizing the results of DeForcrand et al. 
However we have chosen to improve the topological charge via an improved field strength 
tensor. In particular, we employ an 0{a i ) improved definition of F^ v in which the standard 
clover-sum of four lxl Wilson loops lying in the /i, v plane is combined with 2x2 and 



3x3 Wilson loop clovers. Bilson-Thompson et al. [^] find 




h.c. 



(2-4) 

Traceless 



where W nxn is the clover-sum of four nxn Wilson loops and where F^ v is made traceless 
by subtracting 1/3 of the trace from each diagonal element of the 3x3 color matrix. This 
definition reproduces the continuum limit with 0(a 6 ) errors. We employ the plaquette 
measure of the mean link 

u = (~Re irU^fJ^ , (2.5) 

which is updated after each sweep through the lattice. The mean link u rapidly tends to one 
under cooling, thus reproducing the classical limit. However, early in the cooling procedure, 
the fields are not classical and uq serves to tadpole improve the definition of F^ v . 

This improved field-strength tensor can be used directly in Eq. (|1.2j ) resulting in a topo- 
logical charge which is automatically free of discretization errors to the same order as the 
field-strength tensor. On self-dual configurations, this operator produces integer topological 
charge to better than 4 parts in 10 4 . Furthermore, since the gluon (i.e., Yang-Mills) action 
is also based upon the field-strength tensor, it is possible to create what we refer to as a 
reconstructed action based upon the improved field-strength tensor. The value of the action 
calculated with the reconstructed action operator can be compared with the value calculated 
with the standard improved action operator Eq. (|2.3|) at each cooling sweep as a double- 
checking mechanism to ensure that the two different approaches to tree-level improvement 
removing 0(a A ) errors yield consistent results. 

The principal criteria by which we may judge the value of an improvement scheme are 
how quickly and how closely the results for the topological charge Q approach integer values, 
and the stability with which they remain at that integer. On a lattice, the discretization 
errors prevent us from obtaining exactly integer results. However the more successful the 
improvement program, then the more rapidly we can expect cooling to lead us to a stable Q 
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FIG. 1. An example plot of how improved cooling stabilizes the action (circles) and topological 
charge (triangles) at consistent values. The action has been rescaled by dividing it by So, the 
action of a single instanton. 



and the closer this will be to an integer. Figure [l| shows how both the action and topological 
charge approach the same integer value as a function of the number of cooling sweeps. This 
is exactly as we would expect for Q. As we approach self-duality on the 4-torus we appear to 
recover S/Sq = \Q\ just as in the continuum infinite-volume limit. Recall that the positivity 



of {f; v ± 



ensures that < f d 4 x{F" ll/ ± F a f = 2f rf 4 x[(F" ) 2 ± F a F a \ and hence 



S 



hi/ ) 



> 



u djl /us 1 hi/ 



57T 



g z 



Q\=S \Q 



±F a 



(2-6) 



saturate this 



where we have denned S = 8ir / g . Self-dual configurations F^ v 

identity, i.e. , \ Q \ — S/Sq. This result applies independent of the shape of the space-time 
manifold and so applies on the continuum 4-torus as well as for infinite space-time. However, 
note that Nahm's theorem [26| implies that there is no self-dual | Q |= 1 configuration 
possible on the 4-torus. In infinite space-time So is the single instanton or anti-instanton 
action. 

While an in-depth comparison of the various cooling schemes is beyond the scope of 
this report, we summarize by noting that it has been determined [fZIJ that 3-loop improved 
cooling with a 3-loop topological charge operator gives excellent results in terms of how close 
the calculated values of S/S and Q come to integer values, the stability with which they 
remain at that integer, and in terms of the speed with which cooling can be performed. 



III. OVERLAP FERMIONS AND ZERO MODES 



The overlap fermion formalism |27j provides a way of realizing exact chiral symmetry on 
the lattice. The massless overlap-Dirac operator can be written as 
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D o (0) = (l+ 75 e(#w)), 



(3.1) 



where e(H w ) is the matrix sign function 



\H,, 

and D w is the usual Wilson-Dirac Operator 



D w -M-aam,b/3n 



e(H w ) = r^r , # w = 7sAv, (3.2) 



4 

= S atb 8 aj p8 m)n - K - 7M)a/3^( m )a6^m,n-A + (1 + l^a^U^m - fl) ab$m,n+(i\ (3-3) 

(3.4) 

and where k is the hopping parameter, 

k = . (3.5) 

-2m + 8 v ; 

The bare mass parameter for the Wilson kernel in the overlap formalism m has to be in 
the range (0,2) for D o (0) to describe a single massless Dirac fermion. Note that the m used 
in this context is the negative of the bare mass used in simulations with the Wilson quark 
action itself. In principle, any value of m in the above range should give the same continuum 
theory. But on a finite lattice, where volume and lattice spacing are finite, the results for the 
overlap action can depend on the the value chosen for m. For the purposes of using overlap 
fermions to be most sensitive to the topology of the background gauge fields ||12|| , m has to be 
chosen m > mi(g 2 ) for some m\{g 2 ) going to zero in the continuum limit. For typical gauge 
configurations m,i(g 2 ) is slightly less than m c , where m c is the critical value of m at which 
the pion mass extrapolates to zero in a simulation with ordinary Wilson fermions. In the 
hopping parameter formalism, the k should to be in the range (k c , 0.25) as tree level. The 
massless overlap operator D o (0) has been shown to satisfy the Ginsparg- Wilson relation 

{75,^0(0)} = ^0(0)75^0(0). (3.6) 

Its spectral properties are 

• The modulus of the eigenvalue of D o (0) lies in the range [0, 2]. It has exact zero 
eigenvalues which are associated with topology. The corresponding eigenvectors are 
also eigenvectors of 75. They need not occur in pairs. There can be n+ zero-modes 
with eigenvalues of 75 equal to 1 and n_ zeromodes with eigenvalues of 75 equal to -1. 
It also has eigenvalues equal to ±2 which also have definite chirality. Their difference 
{n' + — n'_) is equal to (n_ — n+). 

• Non-zero eigenvalues of D o (0) are complex, their moduli are less than 2, they come as 
pairs and are conjugate to each other. The associated eigenvectors are not eigenvectors 
of 75. 
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In practice, we calculate a small number of low- lying eigenvalues and eigenvectors of 
Dl(0)D o (0). Note that Dl(0)D o (0) commutes with D o (0) and it is hermitian and posi- 
tive definite. It also commutes with 75 and can be simultaneously diagonalized. Hence, 
Dl(0)D o (0) has zero eigenmodes of definite chirality which are also the zero eigenmodes of 
D o (0). 

We compute low-lying eigenmodes of D|(0)Z) o (0) using the Ritz functional algorithm [28 



In the computation of the overlap operator, the time consuming part is the calculation of the 
matrix sign function e(H w ). There are several numerical approaches by which to approximate 
the sign function e(z) |2^-^TJ. We adopt the optimal rational approximation |JIJ with a ratio 
of polynomials of degree 12 in the Remez algorithm. We find the error to the approximation 
of e(z) to be within 10 -6 in the range [0.04, 1.5] of the argument z. To improve the accuracy 
as well as efficiency of computation of the matrix sign function e(H w ), a number of low-lying 
eigenvalues of H w whose absolute value are less than 0.04 are projected out. Since 

Dl(0)D o (0) = Dl(0) + D o (0), Dt(0) = 75^(0)75, (3.7) 

we use chiral states in the Ritz functional algorithm, and hence can save one matrix multi- 
plication |32|] per iteration. 



IV. RESULTS AND COMPARISON 

Configurations have been generated on an 8 3 x 16 lattice at both (3 = 4.80 and (3 = 4.38 
as well as on a 12 3 x 24 lattice at (3 = 4.60. Configurations are selected after iVtherm = 5000 
thermalization sweeps from a cold start and every 7V sep = 500 sweeps thereafter with the 
average link uq fixed at the time of the first sample configuration being taken. Lattice 
parameters are summarized in Table |. 

We first use the gluon field definition to calculate the topological charge Q by using the 
three-loop improved field strength as described in Sec. |1[ For each configuration we 
measure the topological charge by cooling with the 3- loop improved action for just enough 
cooling sweeps that Q is integer to within 1%. This typically requires from 1-30 cooling 
sweeps depending on the lattice spacing and the particular configuration. We retain both 
the original uncooled configurations as well as these "just-cooled" configurations. We denote 
the topological charge obtained from the "just-cooled" gluon field configurations as Q g . 

Secondly, we use the overlap formalism to calculate the index of the gauge configurations 
at two k values, i.e., k = 0.19 corresponding to m — 1.36 and k = 0.2499 corresponding to 
m = 1.999. The second value approaches the largest allowed value of m to describe a single 
massless Dirac fermion. We extract the overlap index for both the original uncooled and the 
"just-cooled" cooled configurations and denote the results by Qd and Qd-cooied respectively. 



TABLE I. Parameters of generated lattices. 



Action 


Volume 


-^therm 


N sep 


P 


a (fm) 


U 


Physical Volume (fm) 


Improved 


8 3 x 16 


5000 


500 


4.80 


0.093 


0.89625 


0.75 3 x 1.50 


Improved 


8 3 x 16 


5000 


500 


4.38 


0.165 


0.87821 


1.32 3 x 2.64 


Improved 


12 3 x 24 


5000 


500 


4.60 


0.125 


0.88888 


1.50 3 x 3.0 
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TABLE II. Results for f3 = 4.80 improved gauge configurations on an 8 3 x 16 lattice with spacing 
a = 0.093 fm: Qd is calculated from the zero-modes of the overlap operator on the original uncooled 
configurations; Q g is obtained using the improved topological charge operator for the "just-cooled" 
configurations, (i.e., configurations after just-enough improved cooling sweeps to bring Q within 
1% of integer); Qd-cool is also calculated from the zero- modes of the overlap operator, but on the 
"just-cooled" configurations. 



1 ' II 

uration # 


Qd{K = 0.19) 


Qd(K = 0.2499) 


Q g 


/ j 

Vd— cool 


1 














2 














3 














4 


+1 


+1 


+1 


+1 


•5 














6 


— 1 


-1 


— 1 


— 1 


7 














8 














9 














10 














11 














12 














13 














14 














15 














16 














17 














18 














19 














20 














21 














22 














23 














24 


-1 


-1 


-1 


-1 


25 


-2 


-2 


-2 


-2 


26 


+1 


+1 








27 














28 














29 














30 
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TABLE III. Results for (3 = 4.38 improved gauge configurations on an 8 3 x 16 lattice with 
spacing a = 0.165 fm: is calculated from the zero- modes of the overlap operator on the original 
uncooled configurations; Q g is obtained using the improved topological charge operator for the 
"just-cooled" configurations; Q g s is obtained using the improved topological charge operator for 
the configurations obtained after just 3 improved cooling sweeps; Qd-cool is also calculated from 
the zero-modes of the overlap operator, but on the "just-cooled" configurations. 



Configuration # 


Q a {k = 0.19) 


Q d (n = 0.2499) 


v B J 




— cool 


1 








+1 








2 


+3 


+3 


+4 


+3 


+3 


3 


+1 














4 

















5 








-1 


-1 


-1 


6 


-1 


-1 


-1 


-1 


-1 


7 


+5 


+4 


+4 


+6 


+6 


8 


-1 


-4 


-2 


-2 


-2 


9 


+1 


+1 


+1 


+1 


+1 


10 








-1 








11 


-2 


-2 


-3 


-3 


-3 


12 


-1 


-1 


-1 


-1 


-1 


13 


+2 


+2 


+3 


+2 


+2 


14 


+3 


-5 


-4 


-4 


-4 


15 


-3 


-3 


-3 


-3 


-3 


16 








-1 








17 





-1 


-1 








18 


-1 


-1 


-1 


-1 


-1 


19 


+1 


+1 


+1 


+1 


+1 


20 


+1 








-2 


-2 



The results for the small- volume 8 3 x 16, (3 = 4.80 lattice are collected in Table |I[ We see 
that since the physical volume of our lattice is rather small, only a few configurations have 
non-trivial topology. The index Q d of the overlap operator calculated at the two different k 
values are the same for all 30 configurations. The range of cooling sweeps needed to get Q to 
within 1% of integer was between 1 and 7 with an average of approximately 3. The index for 
these "just-cooled" configurations differs from these for only one of the thirty configurations, 
i.e., configuration 26. For all of the configurations the gluon topological charge Q g is identical 
to the the index of the overlap operator extracted from "just-cooled" configurations Qd-cooied- 
These results indicate that on such a fine lattice (a = 0.093 fm) topology is relatively well 
represented and the index theorem is "almost" valid for uncooled configurations. We see 
also that for "just-cooled" configurations the index theorem appears to be perfectly satisfied. 
Future investigations should explore the possible volume dependence of these results. 

We next turn to the coarser (a = 0.165 fm), larger- volume 8 3 x 16, (3 — 4.38 lattice. The 
results for this lattice are shown in Table [ITT]. Since the physical volume of the lattice is 
larger there is more nontrivial topology than before, which is reflected in the large Q values. 
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We find that Qd now differs in 30% of cases depending on which value of k is used in the 
overlap kernel. In some cases, e.g., for configuration 14, the disagreement is very significant 
indeed. The range of cooling sweeps needed for this coarse lattice was between 4 and 28 with 
an average of approximately 12.5. In addition to measuring the gluon topological charge 
on these "just-cooled" configurations Q g we have on this lattice also attempted to identify 
the gluon topological charge after just three cooling sweeps, i.e., Q g 3. This was done by 
identifying the nearest integer to the gluon topological charge after just three cooling sweeps 
and was motivated by the observation that the action density and the reconstructed action 
density matched within 10% after three sweeps. We see that Q g3 is reasonably consistent 
with the robust Q g , but there are significant differences. Q g % also has significant differences 
from the Qd values, The disagreement between Qd and Q g are more significant on this coarse 
lattice. 

These differences may be interpreted by considering the size of the topological objects 
giving rise to exact zero modes |]3"S|| . At k — 0.19, the overlap operator will typically to 
miss zero modes associated with small topological objects. Indeed, for the six configurations 
where Q a (k = 0.19) ^ Q a (k = 0.2499), Q d (n = 0.2499) agrees better with Q g3 than Q g . 
This is as one might expect, as further cooling will remove topological objects smaller than 
the dislocation threshold of the improved cooling algorithm, which is typically two lattice 
spacings. 

These results all suggest that for uncooled configurations on this coarser lattice topology 
is not well represented and the index theorem is badly violated. However, from the perfect 
agreement between Q g and Qd-cooi, we see that for "just-cooled" configurations the index 
theorem is again perfectly satisfied. In that sense, we conclude that our "just-cooled" 
configurations are indeed smooth enough. 

We now present in Table [TV] the results for the third lattice with an intermediate lattice 
spacing a = 0.125 fm corresponding to (3 = 4.60 and with lattice size 12 3 x 24. This lattice 
has the largest physical volume. The range of cooling sweeps needed for this medium-spaced 
lattice was between 2 and 14 with an average of approximately 8. As we might anticipate, 
the agreement of Q g with in this case is worse than for the fine lattice but better than 
for the coarse lattice. We do not have results for the large k = 0.2499 case here, since on 
this larger lattice this marginal choice of k proved numerically difficult. Calculating Qd-cooi 
with its "just-cooled" configurations gives perfect agreement with Q g for all configurations. 

It is expected that (Q 2 ) should scale approximately as the volume V for large volumes, 
since the topological susceptibility is given by (Q 2 )/V and should be volume independent. 
The ratio of this 4- volume to that for the coarse lattice is 1.66. The mean Q 2 per configu- 
ration for this lattice is easily seen from Table [TV] to be (Q 2 ) = 7.3, whereas for the coarser, 
smaller- volume lattice in Table |TTT| we find (Q 2 ) = 4.85. We see that (Q 2 )ug/ (Q 2 )smaii = 1-51 
which is approximately equal to Vbi g /V^ ma n = 1.66. This level of agreement seems reasonable 
for such modest volume lattices and small numbers of configurations. 

We also observe that for this largest lattice there appears to be a significant imbalance 
in the sign of Q in the ensemble and in particular the sign is consistently negative for the 
last 8 configurations in the ensemble. If the configurations were uncorrelated, this would be 
very unlikely to occur. This suggests that for larger lattices, when measuring quantities that 
are very sensitive to topology, one should use an increased number of thermalization and 
separation sweeps. To balance the topological charge in the ensemble one should perhaps 
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TABLE IV. Results for (3 = 4.60 improved gauge configurations on an I2 3 x 24 lattice with 
spacing a = 0.125 fm: Qd is calculated from the zero-modes of the overlap operator on the original 
uncooled configurations; Q g is obtained using the improved topological charge operator for the 
"just-cooled" configurations; Qd-cool is also calculated from the zero-modes of the overlap operator, 
but on the "just-cooled" configurations. 



aration # 


Q d {n = 0.19) 


Qg 


Qd— cool 


1 


-1 


-i 


-1 


2 


+1 


+i 


+ 1 


3 


-4 


-4 


-4 


4 











5 


+1 


+1 


+1 


6 


-2 


-2 


-2 


7 


-2 


-2 


-2 


8 


-1 


-1 


-1 


9 











10 


-2 


-2 


-2 


11 


+1 


+1 


+1 


12 











13 


-3 


-3 


-3 


14 


-4 


-3 


-3 


15 


-2 


-4 


-4 


16 


-1 


-2 


-2 


17 


-5 


-5 


-5 


18 


-2 


-3 


-3 


19 


-5 


-5 


-5 


20 


-4 


-4 


-4 
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consider doubling the ensemble size by adding parity-transformed link configurations [34 
which leave the action invariant but reverse the sign of Q. 



V. CONCLUSIONS AND DISCUSSIONS 

We have shown that overlap fermions are suitable for use in the study of topology and 
the Atiyah-Singer index theorem in lattice simulations. We have shown that with an im- 
proved gluon action and an improved definition of the gluon topological charge operator we 
can obtain near integer topological charge (within 1%) within a relatively small number of 
cooling sweeps. The finer the lattice the fewer cooling sweeps are needed to obtain these 
"just-cooled" configurations. We found that for lattice spacings of a = 0.093, 0.125 and 
0.165 fm we needed on average 3, 8, and 12.5 improved cooling sweeps respectively. For all 
configurations on all lattices, we found the index theorem satisfied for these "just-cooled" 
configurations. Even for the finest lattice and with an improved topological charge operator 
some small number of cooling sweeps was needed to obtain an integer gluon topological 
charge. On the finest lattice the index of the overlap operator appeared independent of n 
and agreed with the "just-cooled" overlap index 29 out of 30 times. This lattice appears 
then to be almost fine enough that the index theorem has meaning without cooling. For all 
lattices the index theorem appeared to be fully satisfied for the just-cooled configurations. 
This provides us with a clear benchmark for smoothness and lattice spacing when calculating 
lattice quantities which are very sensitive to topology. 
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